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Abstract 



. It is shown that if the computer model of biological ageing proposed by 

j~^ ■ Stauffer is modified such that the late reproduction is privileged then the 

^D ', Gompertz law of exponential increase of mortality can be retrieved. 
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1 Introduction 

Most computer simulations on biological ageing use the Penna model [l|. This 
model, motivated by the accumulation mutation hypothesis, provides a population 



^ ■ of bit-strings that exhibits many features known for real senescence [0. Recently 

Stauffer [Q in a general review of the Penna model has suggested a simpler alter- 
native with less parameters. This new attempt is based on the postulate of a min- 
imum reproduction age and a maximal genetic lifespan. Only these two numbers 
are transmitted from generation to generation, with certain mutations, by asexual 
reproduction. This simple model shows the basic features required for senescence. 
The age distribution shows an increase of mortality with age. The model reproduces 



successfully the catastrophic senescence of salmon [^ and some specific human as- 
pects like the social structures of technological societies [|] or a transition between 
extinction and survival if a minimal population for survival is required ||^. 

In the following we examine the other modification to the basic model. The 
privilege conditions are introduced to earlier and then to late reproduction. By doing 
this, in case of late reproduction, we achieve the exponential growth of mortality 
in agreement with the Gompertz law. The increase of fertility with age exists in 
nature. Some fish, lobsters and trees share this feature and therefore this life 
strategy has been considered by the Penna model also 0, |^ . 

2 The model 

The population consists of A^ individuals i {i G 1,...,A^) initially. Each individual 
is characterized by three integers: its age a(i), its minimum reproduction age am{i) 
and its maximal genetic lifespan 0^(2) with < am{i) < O'di'i) < 32. The maximal 
lifetime is restricted to 32 time units (called years), the minimum reproduction age 
may be chosen between zero and ad{i) — 1. Within these constraints the values 
of am{i) and a^li) are randomly mutated for an offspring by ±1, away from the 
maternal values, and for each child separately. After an ith individual has reached 
its minimum reproduction age, it gives birth to one offspring with probability b, 
chosen as 

(1) 



Thus during its whole mature life any individual gives one offspring at the cumulated 
probability equal to 1. 

Independently of the genetic death, which happens automatically and unavoid- 
ably if a{i) = ad{i), at each time interval an individual can also die "accidentally", 
with the Verhulst probability N/N^ax- N^ax is called the carrying capacity to ac- 



count for the fact that any given environment can only support populations up to 
some maximal size Nmax- Otherwise the individuals die because of food and space 
limitations. 

In this paper differently than in the standard model we consider populations in 
which the cumulated birth rate is not uniformly distributed among years of repro- 
duction: a^, Om + 1, • • • , fld — 1- Instead we propose to consider: 

(A) privilege to younger third part of a reproduction period: 

jjvoung ^ [ 2^ if am{i) < ^(i) < a^{i) + |Ai ,^. 

\ lb if a^{i) + |Ai < a{i) < ad{i) ^ ' 

(B) privilege to older third part of a reproduction period: 

\ 26 if am{i) + fA^ < a{i) < ad{i) ^ ' 

In the following figures we present characteristics of populations developed when 
three ways of reproduction are considered: standard, (A) and (B) models. To 
avoid possible divergence we follow the Stauffer recipe [^ and shift the birth rate 
accordingly: 

^=^ (4) 

Ai + e ^ ' 

with e = 0.08. In figure 1 we show the population volumes in percentage of the N^ax 
capacity. The population where the reproduction of older individuals is privileged, 
i.e., model (B), is the smallest one and varies in time by 33%. 

Figures 2, 3 compare the distribution of a^ the minimal reproduction age and 
ad the maximal genetic lifespan, respectively, in the three models considered. A sig- 
nificant change in the length of life of individuals is observed when the reproduction 
at old age is favorable, see figure 4 for the distribution of age n{a) in a stationary 
population. 



Finally figure 5 shows the mortality in the populations simulated. The mortality 
is calculated according to the formula g(a) = 1 — n{a + l)/n{a) and plotted on 
log-scale to extract the Gompertz law. 

3 Summary 

In Stauffer's new proposition of a model for biological ageing, by considering the 
privilege conditions for time of reproduction we can easily manipulate the mortality 
in a population. When reproduction at old age is preferred then the model provides 
the mortality in agreement with the Gompertz law. 

In the simple ageing model discussed here the maximal genetic lifespan for each 
individual is limited to 32 years. This limitation was introduced to the model to 
make results better comparable with many Penna model results. But one could ask 
if discarding this assumption would help in getting a nice Gompertz law within the 
standard model. This useful suggestion pointed out to us by the referee has moti- 
vated us for investigations on the simple ageing model. Specially we were interested 
in where the model (B) would drive a population. Preliminary studies show that 
this model simulated without the limitation for the maximal genetic lifespan leads 
to the stationary state such that the distribution of a^ is concentrated around age of 
50. The mortality in this population is presented in Figure 6. One can notice that 
the Gompertz law of the exponential increase of mortality holds for adult individuals 
at age lower than about 35. 
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Figure 1. Usage of environment capacity in time in standard, (A) and (B) models. 

Figure 2. Distribution of minimal reproduction age in standard, (A) and (B) 
models. 

Figure 3. Distribution of lifespan in standard, (A) and (B) models. 

Figure 4. Distribution of age in standard, (A) and (B) models. 

Figure 5. Mortality in populations of standard, (A) and (B) models. 

Figure 6. Mortality in populations of model (B) without the limitation to the 
maximal genetic lifespan. 



0045 - model (A) 




W 




X 

to 




in ,( 



'iiilijii. 



m'^mmimim^'' 



0.040 



c 
o 



5 0.005 



standard 



model (B) 



0.000 J — I — 
2000 



4000 6000 8000 10000 

iteration 



12000 14000 



0.18 



0.16 



0.14 



0.12 



^ 0.10 

CD 
Si 

§ 0.08 



0.06 



0.04 



0.02 



0.00 



distribution of a 



standard 
model (A) 
model (B) 




0.16 



distribution of a^ 



0.14 



0.12 



0.10 



J2 

ro 0.08 

Si 

o 



0.06 



0.04 



0.02 



0.00 



standard 
model (A) 
-▼— model (B) 




0.18 



age distribution 



standard 
model (A) 
model (B) 




age 



0.1 



CO 

t 
o 

E 



0.01 



0.001 




standard 
model (A) 
model (B) 



5 



10 



~\ 

15 

age 



~~r~ 
25 



20 



30 



0.1 



CO 

o 

E 



0.01 



0.001 




age 



